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SUMMARY 


In  conventional  time  domain  beamforming  the  receiver  outputs 
are  delayed  such  that  the  outputs  due  to  a  signal  arriving  from 
a  selected  direction  are  added  in  phase  .  If  the  receiver 
outputs  are  sampled  time  series  then  all  time  delays  must  be 
integer  multiples  of  the  sampling  interval.  In  general  the 
required  delays  are  not  exact  multiples  of  the  sampling 
interval  and  an  appropriate  choice  for  the  approximating  delays 
is  to  choose  the  nearest  integer  multiple  of  the  sampling 
interval. 

This  approximation  of  the  required  time  delays  modifies  the 
polar  response  of  the  array.  In  this  paper  a  model  describing 
the  degradation  in  the  polar  response  of  a  linear  array  of 
equispaced  sensors  at  selected  angles  due  to  the  time  delay 
errors  is  proposed.  At  these  angles  the  perturbed  polar 
response  is  characterized  by  a  number  of  isolated  high 
sidelobes.  These  are  termed  quantization  lobes  and  their 
position  and  approximate  height  are  described  by  the  model. 
Unfortunately  the  model  is  limited  by  its  sensitivity  to 
small  perturbations  in  the  steered  direction. 
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1.  INTRODUCTION 

Time  doaain  beamforming  as  implemented  by  the  DIMUS  (Digital  Multibeam 
Steering)  technique(ref . 1)  involves  the  delaying  of  the  sampled  outputs  of 
individual  receivers.  Whether  these  delays  are  affected  by  shift  register  stages 
or  by  the  use  of  random  access  memories  the  delays  can  only  be  approximated  to 
within  the  nearest  half  sampling  interval.  To  ensure  such  errors  are  kept  to  a 
minimum  the  ratio  of  the  sampling  rate  to  the  highest  frequency  present  should  be 
large.  However,  the  maximum  phase  error  can  still  be  relatively  large,  e.g.  for 
a  frequency  l/8th  of  the  sampling  rate  the  maximum  phase  error  is  */8. 

Anderson (ref . 1)  has  previously  proposed  a  model  which  predicts  the  degradation 
in  the  height  of  main  beam.  His  model  assumes  a  random  distribution  of  the  errors 
of  the  time  delays  with  an  equal  probability  of  lying  within  plus  or  minus  half 
the  sampling  interval  of  the  true  time  delay.  However,  in  many  cases  the  side- 
lobe  response  of  any  array  may  be  an  equally  if  not  more  important  consideration 
than  the  height  or  width  of  the  main  beam  and  is  not  described  by  the  above  model. 

A  model  is  proposed  in  this  paper  whereby  for  certain  steered  angles  the 
perturbed  polar  response  of  a  linear  array  with  equispaced  elements  can  be 
calculated  exactly.  The  model,  described  in  Section  3  is  valid  for  those  angles 
at  which  every  qth  receiver  has  the  exact  time  delay.  Moreover,  for  a  linear 
array  of  N  equispaced  elements  it  is  common  practice  to  steer  N  independent  beams 
which  at  halfwavelength  overlap  at  their  4  dB  points.  If,  for  the  first  beam 
off  broadside  the  progressive  time  delay,  A  r,  is  such  that  the  model  is  applicable 
then  it  can  easily  be  shown  that  the  model  describes  the  polar  response  for  all 
the  other  N  independent  beams. 

A  similar  model  has  been  proposed  by  Allen(ref.2)  to  describe  phase  quantiz¬ 
ation,  and  an  argument  similar  to  that  of  Allen  is  used  in  Section  2  to  relate 
the  mean  square  error  of  a  beam  output  to  that  of  the  time  delays.  The  model 
proposed  by  Allen  considers  only  those  angles  for  which  the  approximating  phase 
delays  have  steps  of  equal  length  and  equal  height.  This  condition  is  relaxed 
in  the  proposed  model  and  as  a  result  a  closer  approximation  to  the  array 
response  is  obtained  for  a  larger  number  of  angles. 

In  Section  4  it  is  shown  how  the  effect  of  the  time  quantization  can  be 
incorporated  as  a  perturbing  factor  on  the  response  of  an  N/q  element  linear 
array. 

The  perturbed  array  response,  at  these  selected  angles,  is  characterised  by 
high  individual  side-lobes.  These  are  termed  quantization  lobes  and  their 
properties  are  discussed  in  detail  in  Section  5.  Finally  in  Section  6  some 
examples  of  the  model  are  given. 


2.  THE  RELATIONSHIP  BETWEEN  TIME  DELAY  AND  PHASE  ERRORS 


Let  Xj(t)  represent  the  time  series  output  of  the  jth  receiver  of  a  linear 

array  of  N  equispaced  receivers.  To  'steer*  a  beam  in  a  direction  0  the  output 
of  the  jth  receiver  should  be  delayed  by  a  r y  where 


_  1 


d  sin  0 


=  ) 


A  r 


d  is  the  separation  of  adjacent  receivers  and  c  is  the  velocity  of  propagation. 

However  if  the  receiver  outputs  are  sampled  in  time  then  it  is  only  possible 
to  delay  the  receiver  outputs  by  integer  multiples  of  the  sampling  interval. 

In  order  to  approximate  the  time  delay  required  for  the  jth  receiver,  an  integer 
m.  is  chosen  such  that 


m.  r0 
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is  niniaised.  It  will  be  shown  in  the  remainder  of  this  section  that  for  N 
large  such  a  choice  minimises  the  mean  square  error  in  the  steered  beam. 
Consider  the  continuous*  beams  y  (t)  and  y  (t)  corresponding  to  exact  and 

C  SL 

approximate  time  delays.  It  follows  that  y  (t)  and  y„(t)  are  given  by: 

C  a 


ye(t) 


N-J. 

=  )  a.  x. (t  -  r .] 

L  3  j  j 


yft)  =  )  a.  x. (t  -  m. 

a  /.  2  r  J 


where  the  a's  are  the  shading  (or  weighting)  coefficients.  Denoting  the  Fourier 
transform  of  (t)  by  (f)  it  follows  that  the  squared  error  of  the  beam  output 

at  time  t,  i.e.,  I  y  (t)  -  y  (t)l 2  is  given  by: 


//  2*  iCf-f'jt 

JJ 


-2wif*  r  . 


Z/2ffifm.T0  2»rifr  .\  / -2irif*  m.r0  -^if'r  \ 

X^XJ^yaja^e  3  -e  3J  (e  ^  -e  J  df  df*  . 


If  the  process  is  stationary  it  follows  that 


<Xj  (f)  Xj^Cf*  )>  =  Sjk  5(f  -  f') 


where <(  )> denotes  the  ensemble  average.  Taking  the  expectation  of  equation  (1), 
substituting  the  above  equation  and  making  the  approximation  of  only  considering 
the  diagonal  terms  in  the  double  summation  it  follows  that  the  mean  square  error 
is  approximated  by: 


f  Ew, 

*  an 


2Jrifm.r<>  2rifr . 


Assuming  identical  receivers  and  a  homogeneous  field  the  above  expression  reduces 
to: 


Eaiaj 


2irifm.ro 

J 


2wifr .  2 


•The  following  argvment  can  also  be  applied  to  sampled  beams  ye(pro)  and  ya(pro) 
provided  a  Fourier  series  expansion  is  used. 


-  3 
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where  S(f)  is  the  power  spectrum  of  the  omnidirectional  field. 
The  expression  (2)  may  further  be  reduced  to: 


2  I  S(f)  E  a?  a^jl  -  cos  2*  f(nuT0  -  T^)i  df 


which  for  2*flm.  t0  -  t.|  small  may  be  reduced  to: 
J  3 


A*1 


Ea*  a.  I  m.T0  -  r  .1 2  [ 

J  J  J  J  / 

J  —  C 


f2  S(f)  df 


Thus  under  the  above  assumptions  the  mean  square  error  of  the  time  delayed  beam 
is  related  to  the  average  error  in  the  approximating  time  delays.  In  particular, 
since  all  quantities  in  the  above  expression  are  positive,  it  follows  that 
choosing  the  m^  such  that  |  a.r0  ~  rj  2S  minimized  is  equivalent  to  minimizing 

the  mean  square  error  of  the  steered  beam. 


3.  FORMULATION  OF  THE  MODEL 

As  shown  in  figure  1  for  a  linear  array  the  exact  time  delays  corresponding 
to  a  particular  steered  direction  are  a  linear  function  of  receiver  number.  For 
selected  directions  the  relative  time  delays  between  certain  receivers  will  be 
integer  multiples  of  the  sampling  rate  and  thus  the  outputs  of  these  receivers 
will  have  the  correct  phase  relationship.  Since  the  array  elements  are  equispaced 
it  also  follows  that  if  the  time  delay  between  the  1st  and  say  the  qth,  receiver 
is  an  exact  multiple  of  the  sampling  interval  then  the  relative  phase  relationship 
between  any  receivers  separated  by  qd  will  also  be  exact.  This  is  illustrated 
by  the  first  example  in  figure  1  where  the  time  delay  required  between  every  second 
receiver  is  equal  to  exactly  one  sampling  interval.  If  the  time  delays  for  the 
odd  receivers  are  chosen  such  that 


I  fflj  t0  -  r  | 

is  minimized  the  stepwise  approximation  to  the  time  delays  as  shown  in  figure  1 
is  obtained.  While  the  approximating  time  delay  between  adjacent  receivers 
is  half  a  sampling  interval  (i.e.  b  To)  in  error  the  relative  delays  of  all  the 
even  numbered  receivers  are  exact  and  similarly  for  the  odd  numbered  receivers. 

In  the  second  example  of  figure  1  the  steered  direction  is  chosen  such  that 
the  relative  time  delay  between  every  fifth  receiver  is  exactly  8  To .  Choosing 
each  time  delay  such  that: 


I  mj  t0  -  r.  | 

is  minimized  it  follows,  similarly  to  the  previous  example,  that  the  relative 
delay  between  any  pair  of  receivers  separated  by  5d  is  also  8  r0 .  This  stepwise 
approximation  to  the  time  delays  is  plotted  in  figure  1.  As  a  consequence  of  the 
above  considerations  the  step  function,  apart  from  a  vertical  shift,  repeats 
itself  with  a  periodicity  of  8  r0 . 
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In  general  there  will  always  exist  angles,  denoted  by  0p  q,  such  that  the 
relative  time  delay  between  any  two  receivers  separated  i>y  a  distance  of  qd  will 
be  an  integer  multiple,  p,  of  the  sampling  interval.  Also,  provided  the  delays 
are  chosen  such  that  I  r0  -  r  |  is  minimized  for  all  j  then  the  stepwise 

approximation  to  the  time  delays,  apart  from  vertical  shifts,  will  be  periodic 
with  a  periodicity  of  p  r0.  This  model  is  a  generalization  of  that  proposed  by 
Allen (ref. 2)  for  phase  quantization.  In  Allen's  model  the  width  of  the  steps 
are  restricted  to  being  equal;  this  restriction  is  removed  in  the  proposed  model 
and  consequently  the  array  response  can  be  evaluated  at  a  larger  number  of  angles. 

At  any  angle  the  progressive  time  delay.  At,  can  be  approximated  arbitrarily 
closely  by  p  T0/q  where  p  and  q  are  integers.  It  then  follows  that  the  delay 
between  every  qth  receiver  will  be  exact  and  due  to  the  periodicity  of  the  time 
delays,  as  discussed  above,  the  array  response  can  be  decomposed  into  a  sum  of 
the  subarrays  with  exact  phase  relationships  between  their  elements.  Before 
discussing  the  effect  of  this  on  the  array  response  some  general  properties  of 
the  factorization  and  time  quantization  will  be  discussed. 

(a)  Finite  Number  of  Receivers 

In  general  the  number,  N,  of  receivers  is  finite  and  so  the  factoriza¬ 
tion  will  only  hold  as  an  approximation  since  q  does  not,  in  general, 
divide  N.  Furthermore  for  finite  N  the  model  will  be  most  useful  when  q 
is  small  compared  with  N.  Also  the  number  of  receivers  with  the  correct 
phase  relationship  (or  time  delay)  to  any  given  receiver  will  be  N/q  and 
for  the  model  to  be  most  successful  the  condition  N/q  >  q  should  hold. 

(b)  Coarseness  of  Quantization 

The  degradation  in  the  array  response  due  to  time  quantization  is  a 
function  of  the  ratio  of  the  frequency  of  the  incident  sine  wave  to  the 
sampling  rate.  The  coarseness  of  quantization,  i?,  is  defined  as  the 
ratio  of  the  f^;  the  frequency  corresponding  to  the  halfwavelength  of  the 
array,  to  l/r0  ;  the  sampling  rate.  Thus 

V  «  fjjTo 

=  c  T0/2d 

where  c  is  the  velocity  of  propagation. 

The  above  intuitive  definition  has  an  alternative  interpretation.  At  a 
frequency  f,  the  beamwidth,  0g,  is  approximated  by 

0B  *  arcsin  (^  f) 

The  delay  between  adjacent  receivers,  A  r,  to  shift  the  main  beam  by  one  beam- 
width  is  then 


At  =  -sinflB 


1_ 

Nf  • 
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Thus,  apart  from  the  1/N  factor,  the  coarseness  of  quantization  is  the  ratio  of 
the  sampling  interval  to  the  delay  between  adjacent  receivers  required  to  shift 
the  main  beam  by  one  beamwidth.  (Note  that  this  can  be  defined  at  any  frequency 
but  for  the  remainder  of  this  paper  only  f^  will  be  considered).  This  is  a 

convenient  measure  of  degradation  as  for  small  tj ,  i.e.  fine  quantization,  the 
approximating  delays  will  be  more  accurate. 


4.  DERIVATION  OF  THE  PERTURBED  ARRAY  RESPONSE 


4.1  Factorization 

Consider  the  array  steered  in  a  direction  such  that  the  relative  delays 
between  receivers  separated  by  a  distance  of  qd  is  an  integer  multiple,  p,  of 
the  sampling  interval.  As  discussed  in  Section  3  all  receivers  separated  by 
multiples  of  qd  will  have  the  correct  phase  relationship.  Consequently  the 
array  response  can  be  factored  into  the  product  of  the  response  of  an  array 
of  N/q  elements  with  the  response  of  an  array  of  q  elements. 

Denote  by  6 the  angle  such  that  the  time  delay  between  adjacent 

receivers  (A  r)  satisfies 


q  A  r  =  p  r0 


qd  sin  0 

_ _ p,q 


c 


(3) 


The  response, 
wave  incident 


P  (0)  when  the  array  is  steered  in  a  direction  6  to  a  sine 
p,q  p,q 

from  a  direction  0  is  given  by 


P 

p.q 


(0) 


J  Y  2»rif(m.  To-(j-l)At)  2 

WL  e 


j=i 


(4) 


where  the  hk  are  chosen  such  that  Iiiu  To  -  t^|  is  minimized  and  A  t  equals 


d  sin  0 


From  the  ’periodicity'  of  the  delays  demonstrated  in  Section  3  it 


follows  that  if  the  delays  for  receivers  1  to  q  are  mi  r0 ,  m2  r0 ,  .  mq  r° 

(where  m^  =  p)  then  the  delays  for  receivers  q  +  1  to  2q  are  (mi  +  p)  t0 , 

Substituting  these  iik's  in  equation  (4)  and 


(nip  +  p)  To  and  so  on. 


summing  over  every  qth  receiver  P^  ^(0)  can  be  written  as: 
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*v  n 

p  (0)  1.  ^  e2jrif  (j-1)  (p  r0-q  A  t)  2k\  f  m,  rc 


+  \  e2’rif^"1^  T°"q  A  ^  e27rif^m2  T°  '  A  t) 


+  a2rif(j-l)  (pro-qit)  e27rif(mq  To'(q"1)At) 


which  reduces  to 


Vq^ 


.  i  -f  .»«0-»)(p  T*> -q  a  t)  f  ■  »“<■*  r'"<k-1’At): 


N*  /  6 


From  equation  (3)  the  above  expression  may  be  further  expressed  as: 

N/q  2ffif (j -1)  3^.  (sin  0  -sin  <p )  q  27rif(m  r0-(k-l)  At) 

pp.,<«  ■  )  *  p-q  )  e 


sin2  Jr  (sin  0  -  sin  <p) 

_ £ _  P  i  q _ 

(N/q)2  sin2fff  —  (sin  6  -  sin  0) 


3  e 


27rif(mk  T0  -  (k-l)At) 


=  “S/q  Rq- 

From  inspection  the  first  factor  R®^  represents  the  exact  polar  diagram  of 

a  linear  array  of  N/q  elements  with  a  distance  qd  between  adjacent  receivers. 

Furthermore  it  is  the  response  of  an  array  steered  in  the  direction  0 

p»q 

The  second  term  R°  is  a  sum  of  q  phase  factors  and  contains  the  effect  of  the 

time  delay  errors  in  the  choice  of  the  ijk  .  R^  will,  in  the  main,  describe 

the  response  of  the  array  and  the  R^  will  be  interpreted  as  a  perturbation 

upon  this  response.  Often  R°  (see  examples)  will  have  a  simple  form  but  in 

general  the  phase  factors  i.e.  the  2rfm^  are  not  linear  in  k  and  as  a 

consequence,  a  simple  expression  for  R°  is  not  possible. 
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In  particular  the  peak  of  k 


will  not  always  lie  in  the  direction  0 

p»q 


as  it 


would  if  the  phases  were  exact.  However,  it  follows  from 
of  Section  2  that  the  distortion  will  be  minimized  since  tl 

most  closely  approximate  to  0^  . 

It  should  also  be  noted  that  as  a  trivial  example  when  p  =  1  the  above 
response  reduces  to  the  exact  polar  response.  This  is  due  to  the  fact  that 
at  certain  selected  angles  each  receiver  delay  is  an  exact  multiple  of  the 
sampling  rate  and  no  distortion  of  the  polar  response  occurs.  This 
illustrates  the  very  strong  dependence  of  the  degradation  on  the  angle  at 
which  beams  are  steered. 

A  more  illuminating  example  is  given  by  considering  the  polar  response  of 
an  array  for  a  steering  angle  0 , ,2  where  alternate  receivers  have  the  correct 
phases.  The  factorization  of  the  array  response  into  the  product  of  the 
response  of  an  N/2  element  array  with  a  two  element  array  can  be  seen 
intuitively  with  reference  to  figure  1.  This  example  is  discussed  in  more 
detail  in  Section  6. 


the  arguments 
e  peak  of  R°  will 

q 


4.2  Discussion  and  interpretation 


The  array  response  at  A/2  will  be  considered  since  the  response  at  any 
other  frequency  can  be  derived  from  it.  Since  the  elements  contributing  to 
R^q  are  separated  by  qd  an  incident  plane  wave  whose  frequency  corresponds 

to  X/2  (=  d)  would  be  spatially  undersampled.  A  consequence  of  this  is  the 
spatial  aliasing  of  the  field  and  the  appearance  of  q  grating  lobes  in  the 
array  response  factor  R^  .  This  is  illustrated  in  figure  2.  To  obtain 

P  (0)  the  modulating  factor  R°  must  be  calculated.  As  discussed  the  peak 

p  j  q  q  q 

response  of  R^  will  most  closely  approximate  that  of  0^  ^  and  a  typical  R 
is  illustrated  in  figure  2.  The  modulating  effect  of 
is  always  less  than  one,  has  two  major  consequences: 


R^  whose  magnitude 


(a)  The  height  of  the  main  peak  will  be  reduced.  However  some  simple 
examples  show  that  unless  the  quantization  is  very  coarse  the  main 
peak  will  be  reduced  by  less  than  a  few  dB. 


(b)  The  height  of  the  grating  lobes  will  be  substantially  reduced.  In 

general,  however,  these  grating  lobes  will  not  be  reduced  to  the  height 
of  the  other  sidelobes  and  will  appear  as  isolated  peaks  of  high  leakage 
into  the  main  beam.  The  number,  position  and  height  of  these  grating 
lobes  can  be  readily  calculated  and  in  the  following  section  some 
expressions  for  these  are  obtained.  To  avoid  confusion  with  the 
grating  lobes  due  to  spatial  under- sampling  when  the  exact  time  delays 
are  used  the  lobes  due  to  the  time  quantization  will  be  termed 
quantization  lobes. 


R° 

q 


is  termed  the  'perturbing  factor'  since  it  contains  the  effect  of  the 


quantization  errors  in  its  non  linear  phase  terms. 


5.  DEGRADATION  IN  ARRAY  PERFORMANCE 


The  degradation  in  array  performance  may  be  quantified  by  use  of  the  model. 
Although  not  every  steered  angle  will  correspond  to,  or  be  approximated  by  a 
0  ,  there  will,  in  general  be  a  sufficiently  large  number  of  steered  angles  such 

p.q 

that  broad  features,  of  the  degradation  may  be  ascertained.  In  particular  the 
dependence  of  features  such  as  the  number,  position  and  height  of  the  quantization 
lobes  and  the  degradation  of  the  main  beam  on  the  parameters  such  as  frequency, 
coarseness  of  quantization  and  angle  of  steer  can  be  determined  by  reference  to 
this  model.  In  the  following  sections  a  quantification  of  these  effects  will  be 
attempted. 
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5.1  Quantization  lobes 

At  an  angle  0^  and  at  a  frequency  corresponding  to  the  half  wavelength 

of  the  array  there  will  be  (q-1)  quantization  lobes.  At  an  arbitrary 

frequency  f  the  number  of  quantization  lobes  is  [  f/f,  (q-1)]  where  [  x]  denotes 

-3 

the  integer  part  of  x. 

The  quantization  lobes  occur  at  angles  0  satisfying: 

in 

(sin  0  -  sin  0  )  =  m 

A.  p,q  m' 


m  1 )  2,  • • • y  q - 1 • 


sin  0  =  -“t  +  sin  0 

m  qd  p,q 


which  at  halfwave length  reduces  to: 


However 


sin  0  =  +  sin  0 

m  q  p,q 


sin  0  =  —  t) 


where  V ,  as  defined  in  Section  3  is  a  measure  of  the  coarseness  of  quanti¬ 
zation.  Consequently  equation  (5)  reduces  to 

sin  0m  =  -PI?  T  for  m  =  1,  2,  ....  q-1.  ( 


0 

The  quantization  lobes  contained  in  R^,^  have  unit  amplitude  and  their  height 

in  the  perturbed  polar  response  is  determined  by  the  amplitude  of  R°  at  the 

o  * 

angles  0  .  Since  R  does  not  in  general  have  an  analytic  form  an  approxim- 

m  o  ^ 

ation  to  R^  is  given  by  the  response  of  a  q  element  array  steered  at  broad¬ 
side.  (See  for  example  figure  1).  This  is  equivalent  to  assuming  the  same 

delay  for  receivers  aq,  aq+1,  . (a+l)q.  Such  an  approximation  is  valid 

for  small  values  of  p,  i.e.  at  angles  steered  around  broadside.  Thus,  in 
decibels  relative  to  unity  the  height  of  the  mth  grating  lobe  is  given  by: 


sin  sin  gn 


v  A 

20  lo8'0  - ^d - 

q  sin  (-7  si 


m  =  1 . 
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At  half-wavelength  sin  6 is  given  by  equation  (6)  and  consequently  the  above 
expression  reduces  to: 


20  log,0 


sin  nprj 


sin?™  -  £) 

'  n  n  ' 


(7) 


An  alternative  approximation  for  R° 


is  to  divide  the  p  sampling  intervals 


equally,  but  to  the  nearest  integer  among  the  q  receivers.  This  results  in 

a  linear  delay  versus  receiver  number  for  the  q  receivers.  (See  also 

figure  1).  This  will  be  a  better  approximation  than  the  previous  one, 
although  for  some  angles  the  two  are  identical.  It  follows  that  R°(0)  is 


approximated  by: 


sin  (sin  6  -  sin  <t>) 
q  sin  ^  (sin  6  -  sin  0) 


where  sin  6  is  given  by 


(8) 


sin  6 


2 


* 


The  height  of  the  mth  quantization  lobe,  hm>  is  given  by  R°(®m),  evaluated 

when  d  =  X/2  and  where  the  assumption  that  this  height  does  not  vary  greatly 
with  frequency  has  been  made. 

Substituting  equation  (6)  for  0m,  the  relative  height  (in  dB)  of  the  mth 
quantization  lobe  h^,  is  given  by: 


h  = 
m 


20  log,0 


sin  w?(q|fcj  -  p) 
q  sin^m^J  -  fj  *  ^ 


(9) 


Note  that  as  the  sampling  rate  increases  then  ij,  the  quantization  measure, 
approaches  zero.  As  i)  •*  o  both  approximations  for  the  heights  of  the 
quantization  lobes  approach  zero  which  would  be  expected. 

5.2  Main  beam 


Figure  2  shows  that  |xcept  for  the  special  case  of  when  the  delays  are 
exact  and  thus  R°  and  Rj^  have  their  main  lobes  pointing  in  the  same 


direction  the  height  of  the  main  beam  in  the  perturbed  array  response  will 
be  reduced.  The  reduction  in  the  height  of  the  main  beam  will  be  a 
function  of  the  wavenumber  separation  of  the  main  beams  of  R°  and  R^  . 


This  in  turn  will  be  determined  by  the  coarseness  of  the  quantization. 

For  example  for  0 j ,2  the  main  beam  of  R°  is  always  directed  towards 

broadside  and  that  of  in  the  desired  steered  direction,  0\  ,2  .  Now  if 


the  quantization  is  very  coarse  then  6 j  ,2  may  correspond  to  endire  directions 
and  consequently,  see  figure  2,  the  reduction  in  the  height  of  the  main  peak 


*  {  xl  has  been  used  to  denote  the  nearest  integer  to  x. 
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would  be  large.  For  fine  to  medium  quantization  6 , ,2  will  correspond  to 
broadside  directions  and  consequently  the  reduction  in  the  height  of  the  main 
peak  will  be  small.  The  reduction  in  the  height  of  the  main  peak,  Ah,  as  a 
function  of  the  coarseness  of  quantization  is  given  by: 


Ah 


20  log,  o 


sin  mj 


2  sin 


2 


To  illustrate  the  dependence  of  Ah  on  the  coarseness  of  quantization  Ah 
has  been  plotted  in  figure  3  for  q  varying  from  0  to  1.  As  tj,  the  coarseness 
of  quantization,  increases  the  reduction  in  the  height  of  the  main  beam 
increases. 

More  generally  equation  (8)  may  be  used  to  approximate  the  perturbing 
factor  and  hence  the  reduction  in  the  height  of  the  main  beam.  Ah,  is 
approximated  by 


Ah  —  1  - 


d  a 

sin  irq  r-  (sin  6  -  sin  0  ) 

a  P»q 


a  sin 


7rd 


A 

Substituting  for  sin  0  and  sin  0 
reduces  to: 


p.q 


and  evaluating  at  halfwavelength  this 


sin  *„({£}  -  E) 

'  Sl"  "(©  -  I)  ' 


In  general  the  delay  approximations  are  such  that 


2-  and  thus 
<1 


Ah  <  1.  For  medium  to  fine  quantization,  rj;  the  ratio  of  the  frequency  to 
the  sampling  rate  will  be  small  and  consequently  for  reasonable  values  of  q 
the  reduction  will  be  less  than  a  few  dB.  Once  again,  for  a  high  sampling 
rate,  »?  -*  0  and  thus  Ah  *  0. 


5.3  Shading 


From  the  previous  sections  it  can  be  seen  that  the  most  serious  source  of 
degradation  in  the  polar  response  of  the  array  is  the  introduction  of 
quantization  lobes.  These  quantization  lobes  are  a  consequence  of  the 
spatial  aliasing  incorporated  in  and  are  unaffected  by  the  array  weights. 


Thus  the  number  and  angular  position  of  these  grating  lobes  are  the  same  for 
a  shaded  array  as  for  an  unshaded  one.  Furthermore  in  most  schemes  of 
amplitude  weighting  the  shading  coefficients  vary  smoothly  from  receiver  to 
receiver  and  thus  may  bg  considered  to  be  approximately  constant  over  the 
subarray  described  by  R  .  It  thus  follows  that  R°  for  a  shaded  array  is 

^  o  ^ 

approximately  equal  to  the  R^  for  an  unshaded  array.  Thus  the  height  of  the 


quantization  lobes  are  to  a  very  good  approximation  independent  of  the  array 
shading  scheme  used. 

As  would  be  expected  other  factors  such  as  the  beamwidth  of  the  quantization 
lobes  and  the  height  of  the  subsidiary  sidelobes  are  determined  by  the  array 
shading. 
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6.  EXAMPLES 

Consider  a  linear  array  of  N  equispaced  elements  where  f,  is  the  frequency 

corresponding  to  halfwavelength  spacing  and  f0  is  the  sampling  interval.  If  d 
is  the  separation  of  adjacent  receivers  then  the  angles  0  as  defined  in 

Pi*l 

Section  4  satisfy: 


sin  0 


p.q 


pc  To 


where  c  is  the  velocity  of  propagation.  As  an  illustration  of  the  model  the 
following  angles  are  considered. 

6.1  0,,j 


At  this  angle  the  time  quantization  is  such  that  the  delay  between 
adjacent  receivers  is  exactly  one  half  the  sampling  interval.  Consequently 
all  even  and  all  odd  receivers  will  have  the  correct  phase  relationship 
between  them.  The  time  delays  are  shown  in  figure  1  as  a  function  of 
receiver  number.  The  polar  response.  Pi ,2(0),  is  given  by: 


Pi  .2»)  - 


*Nd  a  .  .. 

sin  -y  (sm  0|  ,2  -  sin  0) 

N  2xd  ,  .  a 
^  j  sin  (sin  0 1 ,2  -  sm  4 


2  MTd  •  ^ 

cos  I- ^  sin  0 


The  first  term  is  the  exact  polar  response  of  an  array  of  N/2  equispaced 
receivers*  with  adjacent  receivers  separated  by  2d.  The  N/2  element  array 
is  steered  in  the  direction  0i  ,2 .  The  second  factor,  R? , 


sin  0 


describes  the  polar  response  of  a  dipole.  The  effect  of  the  time 
quantization  is  to  steer  this  dipole  in  the  broadside  direction  for  all  values 
of  0i,2.  If  the  sampling  rate  is  sufficiently  high,  i.e.  the  quantization 
is  fine,  then  0i ,2  will  necessarily  correspond  to  a  direction  near  broadside 
and  little  degradation  will  occur.  However  for  a  low  sampling  rate  0t ,2  will 
approach  endfire  and  so  considerable  distortion  could  be  expected. 

At  f,  the  spacing  of  the  N/2  element  array  is  X  and  so  an  aliased  lobe  will 

appear.  Using  equation  (6)  the  position  of  this  quantization  lobe,  in  degrees 
relative  to  broadside,  is  given  by 


sin  0i  =  17  -  1 . 

The  height  of  this  lobe  is  given  by: 


20  logio 


sin 


*3L 

2 


which  is  a  special  case  of  equation  (7) . 


*  For  N  odd  the  above  approximation  will  still  be  valid  provided  N  is 

BtfficifPtly  large. _  J 
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The  perturbed  response  corresponding  to  N  *  30,  =  320  Hz  and  a  sampling 

rate  of  2048  Hz  is  plotted  in  figure  4.  The  polar  diagrams  in  this,  and  all 
subsequent  figures  were  calculated  by  choosing  the  delay  for  the  jth  receiver 
so  as  to  minimize 


nij  T0 


as  discussed  in  Section  2.  The  coarseness  of  quantization,  q ,  has  the  value 
5/32.  o 

The  angle  0j ,2  is  8.9  ,  and  from  figure  4  the  position  of  the  quantization 
lobe  is  -56  and  its  height  is  -12.4  dB  relative  to  the  peak.  Using 
equations  (6  )  and  (9  )  the  predicted  position  and  height  of  the  lobe  are 
-58  and  -12.3  dB  respectively,  and  show  good  agreement  with  the  results  of 
figure  4. 

6.2  0g  ,s 

This  angle  corresponds  to  every  fifth  receiver  having  the  correct  time 
delays  relative  to  each  other.  The  delays  for  each  subset  of  five  receivers 
are  chosen  so  as  to  most  closely  approximate  the  linear  delay  relationship 
and  are  also  shown  in  figure  1.  Assuming  the  total  number  of  receivers  to 
be  a  multiple  of  five  the  perturbed  array  response  is  given  by: 


P8.sC0)  = 


Nfd  ,  .  a  .  j 

sin  -j—  (sin  08  ,s  -  sin  0) 


5-  sin  it  ^  (sin  08  ,s  -  sin  $) 


Rs°(*) 


where  from  figure  1,  R5 (0)  is  given  by: 


R?(*)  =  §5 


1  z  +  zJ  e ^  *  z3  e1^  +  z4  e2i^j * 


where 


2»if (2r0  - 


z  =  e 


and  -  2»ifr<> .  Once  again  the  effects  of  the  time  quantization  are 
incorporated  in  the  R°  factor.  At  a  frequency  of  f,  the  spacing  of  the  N/5 

q  *3 

5X 

element  array  is  and  consequently  four  quantization  lobes  will  appear. 

These,  in  contrast  to  the  exact  array  response  do  not  appear  at  angles 
corresponding  to  the  nulls  of  R?(0}  and  consequently  have  a  finite  height. 

In  figure  5  the  polar  diagram  for  N  =  30,  f^  =  320  Hz  and  r0  =  1/2048  s  is 

plotted  where  0s ,5  is  30°.  The  position  and  heights  of  the  quantization 
lobes  appearing  in  this  figure  are  tabulated  in  the  table  below.  In  order 
to  calculate  the  predicted  values  of  these  lobes,  the  Rg(^)  term,  as 
suggested  in  Section  5,  has  been  approximated  by: 


13  - 
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R?  =  +  *  *z'*W 


f!Sd 


sin  ry-  (sin  0  -  sin 


25  sin*  (sin  d  -  sin  6) 


hi 


y 


where  sin  0  =  4*1 .  Using  equation  (6)  the  position  of  the  mth  quantization 
lobe  is  given  by: 


sin  0. 


16*1  2m 

1"  '  T 


*  j  -  for  m  =  -1,  1,  2  and  3. 


Equation  (9)  for  the  relative  height  of  the  mth  quantization  lobe  reduces  to 


h  =  20  logio 

m 


sin  »2rj 


5  sin  ff 


2n»m 

5 


In  the  table  these  values  have  been  calculated  from  the  above  formulae  and 
are  compared  with  those  measured  from  figure  5.  The  angular  positions  of 
the  quantization  lobes  are  in  good  agreement  whilst  their  heights  only  show 
agreement  at  angles  well  removed  from  the  main  beam. 

The  position  and  heights  of  the  quantization  lobes  for  0j6,s  are  also 
tabulated  in  the  table.  The  polar  response  for  0|6,s  is  shown  in  figure  0 
where  the  same  parameters  as  for  figures  4  and  5  have  been  used.  Once  again 
the  predicted  and  observed  quantization  lobe  parameters  show  good  agreement. 


TABLE  1.  POSITION  AND  HEIGHT  OF  QUANTIZATION  LOBES 


m 

Measured 

Predicted 

Position 

(degrees) 

Height* 

Position 

(degrees) 

Height* 

0a  ,s 

-1 

66 

-19 

64 

-8 

1 

7 

-17.5 

6 

-13 

2 

-17 

-17 

-17 

-15.5 

3 

-43 

-14.5 

-44 

-14.5 

016  >5 

-1 

36 

-14 

37 

-15 

1 

12 

-19 

11.5 

-20 

2 

-10 

-20 

-11.5 

-20 

3 

-37 

-16 

-37 

-17 

*Height  in  decibels  relative  to  unity 
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6.3  Shading  and  angular  perturbations 

In  figure  7  the  perturbed  array  response  corresponding  to  figure  5  but 
with  a  Hanning  window  incorporated  is  plotted.  The  fact,  discussed  in 
Section  5,  that  shading  only  significantly  affects  the  width  of  the  quantiz¬ 
ation  lobes  but  not  their  nuaber,  angular  position  or  heights  is  well 
illustrated  by  this  figure. 

Unfortunately  the  values  of  p  and  q  may  radically  change  for  a  small 
perturbation  of  the  steered  angle.  In  the  above  example  if  the  steered 
angle  is  varied  by  2°  to  28°  then  the  values  of  p  and  q  are  altered  to  3  and 
2  respectively.  As  illustrated  by  figure  8,  the  array  response  correspond¬ 
ing  to  this  new  angle  will  have  drastically  different  characteristics. 

Also  the  examples  given  have  been  chosen  to  show  well  defined  quantization 
lobes.  However  not  all  angles  will  correspond  to  reasonable  values  of  p 
and  q.  The  polar  response  of  the  line  array  of  the  previous  examples 
is  plotted  in  figure  9  for  an  arbitrary  angle  of  34°.  The  grating  lobes  are 
no  longer  so  well  defined  and  the  response  is  characterized  by  a  general 
increase  in  the  sidelobe  level.  This  is  perhaps  the  most  serious  short¬ 
coming  of  the  model. 


7.  CONCLUSIONS 

A  model  has  been  proposed  to  predict  the  degradation  in  the  response  of  an  N 
element  linear  array  at  selected  angles  due  to  the  finite  sampling  interval. 

It  is  capable  of  predicting  the  number,  the  angular  position  and  approximate 
height  of  a  number  of  high  isolated  side-lobes  introduced  by  the  time  delay  errors. 
Furthermore  the  dependence  of  these  quantization  lobes  on  the  parameters  of 
frequency,  sampling  rate  and  array  shading  can  be  deduced  from  the  model. 

Unfortunately  the  model  is  sensitive  to  relatively  small  perturbations  in  the 
angle  of  steer. 
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Figure  5.  Perturbed  response  at  30°  (08 ,s ) 
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Figure  6.  Perturbed  response  at  90°  (016,s) 
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Figure  8.  Perturbed  response  at  28°  (03 ,2 ) 
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